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Abstract 

Using the density-functional perturbation theory with structural optimization, we investigate 
the electronic structure, phonon spectra, and superconductivity of BiS2-based layered compounds 
LaOi_j:Fa;BiS2. For LaOo.sFo. 56182, the calculated electron-phonon coupling constant is equal to 
A = 0.8, and obtained Tc — 9.1 K is very close to its experimental value, indicating that it is a 
conventional electron-phonon superconductor. 

PACS numbers: 74.20.Pq, 74.25.Kc, 63.20.kd 
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The recent discovery of superconductivity in the BiS2-based compounds has attracted 
much attention due to their layered crystal structure and physical properties similar to 
cuprate and Fe-based superconductors. The 6140483 superconductor was found to exhibit 
metallic transport behavior and to show a zero-resistivity state below 4.5 K.^ Soon after, an- 
other type of BiS2-based systems, ReOi_a;Fa;BiS2 (Re = La, Ce, Pr, and Nd), were reported 
to be superconducting with transition temperatures Tc equal to 10.6, 3.0, 5.5, and 5.6 K, 
respectively.-"- For 3140483 (ReOi_a;Fa;BiS2), superconductivity can be obtained by elec- 
trons doping into the insulating parent compound BigOgSs (ReOBiS2). Experimental and 
theoretical studies indicated that these materials exhibit multiband behavior with dominant 
electron carriers originating mainly from the Bi Qp orbitals.-*^"— 

In this Letter we report the calculated results for the electronic structure of LaOi_2:F2;BiS2 
(x = 0, 0.5, and 1), phonon dispersions and electron-phonon coupling of superconducting 
LaOo.5Fo.5BiS2. From the first-principles calculation, it is shown that the parent compound 
LaOBiS2 is a band insulator with an energy gap of ~ 0.8 eV, similar to the ground state of 
cuprate and pnictide superconductors. Upon F doping, the Fermi level moves up and the 
system appears a metallic nature. We obtain electron-phonon coupling constant A = 0.8 and 
superconducting transition temperature Tc ~ 9.1 K in LaOo.5Fo.5BiS2, indicating that this 
compound can be explained as a conventional phonon-mediated superconductor, quite differ- 
ent from the Fe-based materials^*^ but somewhat similar to the Ni-based compounds.— 

We used the full-potential linearized augmented plane wave method implemented in the 
WIEN2K package^^. for the lattice parameter optimization and electronic structure calcula- 
tions. The generalized gradient approximation (GGA) presented by Wu and Cohenr^ was 
used for the exchange-correlation energy calculations. This GGA is a nonempirical approx- 
imation that gives significant improvements of calculations for lattice constants and crystal 
structures. The phonon dispersion and the electron-phonon coupling are calculated using 
density functional perturbation theory (DFPT). The phonon dispersions were obtained in 
linear response via the QUANTUM-ESPRESSO program'^ and ultrasoft pseudopotentials. 
The GGA of Perdew-Burke-Ernzerhof^° was adopted for the exchange-correlation potentials 
with a cutoff of 60 Ry for the wave functions and 600 Ry for the charge density. 

LaOi_a;Fa;BiS2 with x=0 or 1 is crystallized in the ZrCuSiAs type tetragonal structure 
(P4/nmm) with La, Bi, SI in the Bi planes, S2 in between Bi-Sl and La-O(F) planes, and 
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TABLE I: Various optimized structural parameters for LaOi_2,FaBiS2 with x=0, 0.5, and 1. 





LaOBiS2 


LaOo.5Fo.5BiS2 


LaFBiS2 


a (A) 


4.0677 


4.1091 


4.1524 


c(A) 


14.3085 


13.4196 


12.9554 


V (A3) 


236.75 


226.59 


223.38 


ZLa 


0.0890 


0.1073 


0.1276 


ZBi 


0.6309 


0.6145 


0.6090 


zsi 


0.3945 


0.3844 


0.3691 


ZS2 


0.8073 


0.8128 


0.8163 




2.8993 


2.9056 


2.9499 


dBi_52(A) 


2.5231 


2.6621 


2.6862 



0(F) being at the positions 2c(0.5, 0, ZLa), 2c(0.5, 0, zst), 2c(0.5, 0, zsi), 2c(0.5, 0, 252), 
and 2a(0, 0, 0), respectively. The F substitution at the O sites in LaOo.5Fo.5BiS2 reduces 
the symmetry to a space group of P-4m2. Starting from the lattice parameters (a = 4.0527 
A and c = 13.3237 A) reported by Mizuguchi et al.^, we have performed the full structural 
optimization including the lattice parameters and atomic positions, the obtained results 
being summarized in Table [B It is found that, with doping F from a; = to x = 1, the a 
value increases gradually but the c value decreases, accompanied with a decrease of primitive 
cell volume V. By making a comparison in the atomic bond length between LaOo.5Fo.5BiS2 
and LaOBiS2, one finds that the Bi-S2 bond expands due to F doping while dsi-si for the 
Bi-Sl bond remains almost unchanged. 

The calculated band structures and electronic density of states (DOS) are shown in Fig. 
[H The contribution of orbital states to the band structure are characterized by different 
colors: blue (Bi-p), red (Sl-p), and green (0-p and S2-p). For LaOBiS2, the Fermi level 
is located at the upper edge of valence bands with an energy gap of ~ 0.8 eV, indicating 
an insulating behavior. The valance bands spread from around —6 eV to and consist of 
the p states of the O and S atoms. While the La-c? and La-/ states lie above 4 eV, far 
away from the Fermi level, the Bi-p and Sl-p states dominate the conduction bands. For 
doped LaOo.5Fo.5BiS2, the Fermi level enters the conduction bands coming from the Bi-p 
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TABLE II: Phonon mode frequencies (cm ^) at F and M points in LaOo.5Fo.5BiS2. I: infrared 
active, R: Raman active. 



r(o, 0, 0) 


B2{I + R) 





65.4 


145.4 


249.4 270.1 412.2 


E{I + R) 





39.3 


44.0 


78.4 122.8 




153.6 171.3 


184.1 


243.8 345.1 


Ai{R) 


74.8 


138.3 


168.4 


305.9 


M(0.5, 0.5, 0) 




61.3 


81.3 


155.8 


195.8 




55.5 


101.7 


123.5 


189.3 275.8 


Bi 


60.8 


81.6 


169.5 


199.6 


B2 


56.3 


94.3 


103.3 


191.0 299.6 


E 


71.0 


74.9 


148.4 


220.5 269.8 403.4 



and Sl-p states, so that the Bi-Sl layers dominantly contribute to the electronic conduction. 
As shown in Fig. [T|d, the present Fermi level crosses four bands of electron pockets, and 
so LaOo.5Fo.5BiS2 has electron carriers, which is well consistent with the Hall resistance 
measurements.— The DOS at the Fermi level is equal to N{Ep) = 1.22 eV~^ on a per 
formula unit both spins basis, and the corresponding bare susceptibility and specific heat 
coefficient are given by Xo = 4.0 x 10"^ emu/mol and 70 = 3.0 mJ/mol K^, respectively. For 
LaFBiS2, the Fermi level further enters the conduction bands and the energy gap increases 
to ~ 1.1 eV. Meanwhile, the DOS at the Fermi level increases to N{Ef) = 1.84 eV~^ per 
formula unit both spins, yielding xo = 6.0 x 10~^ emu/mol and 70 = 4.5 mJ/mol K^. 

The calculated phonon dispersions of LaOo.5Fo.5BiS2 are plotted in Fig. [2^, in which 
there are 30 phonon bands extending up to ~ 450 cm"^ and the point group at the F and 
M points is D2d- The F and M-point modes can be decomposed as F = 6B2 lO-E 04Ai 
and M = 45i 4Ai 5^2 5^2 6E. Note that the E mode is degenerate, and its 
frequencies are listed in Table [Til The Raman modes at the zone center can be measured 
directly by experiments in near future. The corresponding phonon density of states (PDOS), 
Eliashberg spectral function a'^F^u), electron-phonon coupling A(w) and projected PDOS 
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are shown in Fig. |3l There are three distinct peaks centered at ~ 70 cm~^ {Ai mode at F 
point), ~ 190 cm~^ {E mode at F point), and ~ 250 cm^^ {B2 mode at F point), respectively, 
corresponding to three thick hnes from lower to upper in Fig. [2^. As shown in Fig. the Ai 
mode corresponds to the vertical vibrations with the upper and lower BiS2-layers (La-layers) 
beating against each other. The E mode corresponds to the in-plane stretching vibrations 
with motion along the x{y) axis. In the B2 mode, the La, Bi, and SI atoms vibrate along the 
same z direction, opposite to the vibration direction of the S2, O, and F atoms. From Fig. 
[3l one can see that the bands around 70 cm^^ comes mainly from the PDOS contributions 
of the Bi and La motion, the peak at 190 cm~^ mainly corresponds to the S character, and 
that at 250 cm^^ corresponds to F, O and S character. The calculated Eliashberg spectral 
function has peaks in consistent with the PDOS at the low-frequency regime, and is evidently 
enhanced in the intermediate frequency regime. In the DFPT calculations, the Eliashberg 
spectral function depends directly on the electron-phonon matrix element: 



X(5(£:„k)5(£:m(k+q))'^(w - ^^qu)- 
Here, is the number of k points used in the summation, N{Ep) is the density of states at 
the Fermi level, and Ugi, are the phonon frequencies. The electron-phonon matrix element | 



the electron-phonon interaction of Bi-Sl bond plays an important role in the enhancement of 
a^-F(w). From the A((X') curves in Fig. [3^, where \{uj) = 2 f^[a'^E{Q)/Q]dQ, one finds that 
A(220) ^ 0.7 has been close to A(oo) ~ 0.8, indicating that the phonon modes in the low and 
intermediate frequency regimes have the main contribution to the electron-phonon coupling. 
From the a'^F{uj) spectra in Fig. [3], we get the electron-phonon coupling A = 0.8 and the 
logarithmically averaged frequency uin =192 K. Using further the Allen-Dynes formula with 
the Coulomb parameter fi*= 0.1, we finally obtain = 9.1 K, which is found very close to 
the experimental value Tc = 10.6 K.^ It then follows that the doped LaOo.5Fo.5BiS2 is very 
likely a conventional electron-phonon superconductor. 

In summary, the present first-principles calculations indicate that LaOi_2^^F2:BiS2 changes 
from a band insulator at a; = to a band metal upon doping. Similar to LaFeAsOi-^^Fa;, 
electrons in LaOi_a;F3.BiS2 are doped from the LaO(F) redox layer to the superconducting 
BiS2 layers. The p-orbital states of Bi-Sl layers play an important role in the electron 




(1) 



^nk,m(k+q) 



^ is defined by the variation in the self-consistent crystal potential. It follows that 
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transformation due to their domination of the DOS at the Fermi level. From the phonon 
calculations for LaOo.5Fo.5BiS2, we obtain a strong electron-phonon coupling of A = 0.8 and 
a superconducting transition temperature of = 9.1 K, the latter being very close to its 
experimental value. It then follows that LaOo.sFo. 56182 is a conventional electron-phonon 
superconductor. 

For LaOi_a;Fa;BiS2, we have shown that the electronic states near the Fermi level are dom- 
inantly from the hybridized p orbital states of the Bi and SI atoms, and the electron-phonon 
coupling comes mainly from the contribution of the Bi and S phonons. These facts indicate 
that the BiS2 layers play an important role in the transport and superconducting properties, 
similar to the CUO2 layers of the cuprates or Fe-As layers of the Fe-based superconductors. 
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Figure Caption 



Figure 1. (Color online) Calculated band structure and density of states of LaOi_-rF^BiS2 
for x—0 (a), 0.5 (b), and 1 (c), with orbital character indicated by different colors: Bi-p 
(blue), Sl-p (red), and 0-p and S2-p (green). 

Figure 2. (Color online) (a)CalcuIatcd phonon dispersion curves of LaOo.5Fo.5BiS2. 
(b)Atom vibration pattern of the Ai, E and B2 modes at F point. 

Figure 3. (Color onhne) (a) Phonon density of states G(a;), electron- phonon spectral 
function a^F^u) and electron-phonon coupling A(a;) for LaOo.5Fo.5BiS2. (b) Projected 
phonon density of states for each atom. 
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FIG. 1: Calculated band structure and density of states for LaOi_a;Fa;.BiS2 with x=0 (a), 0.5 (b), 
and 1 (c) with orbital character indicated by different colors: Bi-p (blue), Sl-p (red), and 0-p and 
S2-p (green). 
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FIG. 2: (a) Calculated phonon dispersion curves of LaOo.5Fo.5BiS2- (b) Atom vibration pattern 
of the Ai, E and B2 modes at F point. 
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